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The layer interdependence of transport in an undoped electron-hole bilayer (uEHBL) device was 
studied as a function of carrier density, interlayer electric field, and temperature. The uEHBL device 
consisted of a density tunable, independently contacted two-dimensional electron gas (2DEG) and 
two-dimensional hole gas (2DHG) induced via field effect in distinct GaAs quantum wells separated 
by a 30 nm Alo.9Gao.1As barrier. Transport measurements were made simultaneously on each layer 
using the van der Pauw method. An increase in 2DHG mobility with increasing 2DEG density was 
, observed, while the 2DEG mobility showed negligible dependence on the 2DHG density. Decreasing 

the interlayer electric-field and thereby increasing interlayer separation also increased the 2DHG 
mobility with negligible effects on the 2DEG mobility. The change in interlayer separation as 
interlayer electric-field changed was estimated using 2DHG Coulomb drag measurements. The 
results were consistent with mobility of each layer being only indirectly dependent on the adjacent 
layer density and dominated by background impurity scattering. Temperature dependencies were 
also determined for the resistivity of each layer. 
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03 ! I. INTRODUCTION 

Interest in electron-hole bilayers necessarily arose from the prospect of observing Bose-Einstein condenstion (BEC) 
of excitons in semiconductor double quantum well systems 1 -^ and significant progress towards this goal has been 
made. 3,4 This trend in bilayer research centered on the behavior of the electron-hole pair. The work presented in this 
paper, however, focuses on the individual transport in each layer and the layer interdependence. The latter is the 
primary question we seek to answer here and for which a bilayer device is singularly, exceptionally suited; to what 
extent will the mere presence of a nearby 2DEG affect the transport in a 2DHG and vice-versa? 

The general transport properties of the 2DEG (or 2DHG) system in modulation-doped heterostructures were well 
established decades ago££ Exploiting Coulomb scattering's dependence on the shape of the wavefunction, Hirakawa 
et al. and others demonstrated later that the 2DEG's mobility dependence on density could be altered by deforming 
the wavefunction using external fields from gates Calculations by Kurobe showed that the 2DEG wavefunction 
can be squeezed by making the back-gate voltage more negative, which narrows the 2DEG's confining potential 
I , and moves the wavefunction towards the interface. 10 Furthermore, he demonstrated that squeezing the wavefunction 
■ reduced remote and space impurities scattering times, but enhanced the channel impurity scattering time. To eliminate 
£N| \ remote-impurity scattering, which typically limits the mobility in modulation-doped heterojunctions, and more closely 
inspect the roles of other scattering mechanisms, similar studies were also done using undoped, inverted semiconductor- 
insulator-semiconductor or ISIS structures, pioneered by Meirav et al^ From these studies it was determined that 
background channel impurities dominate scattering at low densities, while interface roughness dominates at higher 
densities j 12 ! 13 More recently, Das Sarma et al. and others showed that background impurity scattering in GaAs 
heterostructures is rudimentary to the 2D metal-insulator transition, which occurs as density is reduced and screening 
of the random potential landscape, due to these impurities, becomes progressively weaken 14 ' 15 These studies were all 
on unipolar devices and, to the author's knowledge, an experimental study of whether general 2D transport is affected 
by a 2D system of opposite charge in close proximity has not been reported. 
' In this paper, the results of an investigation into the layer interdependence of transport in a uEHBL are presented. 
The uEHBL device under study consists of a density tunable, independently contacted 2DEG and 2DHG induced via 
field effect in distinct GaAs quantum wells separated by a 30 nm Alo.gGao.i As barrier^ To populate the undoped wells 
an interlayer electric field E/^ is necessarily established to account for the energy difference between the conduction 
and valence bands. This design affords the following advantages: (1) independent contacts allow for simultaneous 
transport measurements of each layer and Coulomb drag measurements; (2) a tunable density 2DEG and 2DHG 
allows for these measurements to be made as functions of the density in each layer, n and p; (3) an undoped structure 
reduces or eliminates scattering by remote ionized impurities; and, finally, (4) for the same densities, the distance 
between the 2DEG and 2DHG or interlayer separation d can be varied by changing En, and both gate voltages. 

The investigation included mobility and resistivity measurements measured in each layer as functions of n and 
p, E/i, and temperature T. The results indicated that 2DHG transport changed by varying n or En,, while the 
2DEG was largely immune to similar changes in p or En,. Coulomb drag measurements were used to estimate the 
change in d as En, was varied. These results demonstrate that transport in an uEHBL has an asymmetric layer 
interdependence, with only the 2DHG transport having an apparent dependence on the adjacent layer density, n. 
However, this effect appears to be only indirect since increasing n necessitated an increase in the tilting of the 2DHG 
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FIG. 1: (a) Schematic of the uEHBL device, (b) Cartoon schematic of the energy level diagram showing relative locations of 
Vtg, V/l, and Vbg 



confinement potential as well. Presumably, this tilting squeezed the hole wavefunction, and thereby reduced the 
background impurity scattering. The temperature sweeps of the resistivity of each layer also reflected qualitatively 
unique behaviors, with the 2DHG displaying stronger changes at similar densities then the 2DEG. 



II. MATERIAL AND FABRICATION 



A full account of the design and fabrication of the device used in this study was given ini 6 -. The uEHBL device was 
formed from molecular beam epitaxially (MBE) grown GaAs/AlGaAs double quantum well material (wafer EA1286). 
A side profile of the device after full processing is shown in Fig. QJ,. The top and bottom 18 nm GaAs quantum wells 
were separated by a 30 nm Alo.9Gao.1As barrier. Above the top quantum well is a 200 nm Alo.3Gao.7As cladding layer 
and a 60 nm n+ GaAs cap layer. Beneath the bottom quantum well is a 125 nm Alo.3Gao.7As cladding layer and a 
310 nm growth superlattice. Beneath that is a 15 nm GaAs layer, which acts as the second etch stop during backside 
processing and effectively becomes the cap layer. As this layer gets exposed to air, the Fermi level at its surface is 
expected to be pinned at mid-gap. To process the uEHBL device a ~ 25 cm 2 piece of this material was cleaved from 
the wafer and mesa-etched in the shape of a Hall bar with 5 arms extending from each side. 

The carriers in each quantum well were induced via external fields generated by gates on the top and bottom of 
the device, as depicted in[IJa). Each gate covered the central Hall bar region and extended along the entire length 
of half the arms in a geometry which allows for both H xx and Hall measurements on both wells, independently. The 
gate contact pads were positioned at opposite ends of the long axis of the central Hall bar region. 

The only intentional dopants present in the material were in the n+ cap layer, which was also used to form the 
top gate. These dopants do not populate the upper quantum well. A shallow annealing of PdGeAu metal was used 
to make ohmic contact to the top gate. Directly adjacent to the edge of the top-gate at each arm's end were self- 
alignecUi, n-type NiGeAu ohmic contacts. These contacts served a dual role as reservoirs supplying electrons to the 
top quantum well and as the n-type ohmic contacts for transport measurements on the 2DEG. At the ends of the 
remaining arms, AuBe p-type ohmic contacts were formed. 

To process the bottom gate the device underwent the epoxy bond and stop-etch (EBASE) technique 18 . This 
technique entails epoxying the sample to a host substrate topside down and removing material (the original substrate 
and etch-stop layer) down to the growth superlattice by lapping and etching. Once EBASE finished, SiN was deposited 
over the entire mesa surface by plasma enhanced chemical vapor deposition (PECVD) and vias were formed through 
the epi-material to contact the ohmic pads, now adjacent to the host substrate material. 

The back gate metal, TiAu, was then deposited over-top the SiN. The back gate covers the central Hall bar region, 
the five mesa arms attached to p-type ohmic contacts and a small area of these contacts. In this so-called overlap 
configuration, holes are pulled by the backgate into the bottom well from the p-type ohmic contacts^; this is analogous 
to the n-type ohmic contacts' role described above. While both quantum wells are physically in contact with the 
n-type and p-type metal contacts, the gates and mesa configuration is such that only one type carrier is induced in 
each well. 
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III. EXPERIMENT 

A schematic of the energy band diagram of the uEHBL during typical operation is given in Fig. [TJd. To simultane- 
ously establish a 2DEG and 2DHG in the uEHBL devices three different, negative bias voltages, top-gate bias Vtg, 
bottom-gate bias Vbg and interlayer bias Vil, are necessarily used. All these voltages are referenced to ground and 
at least one 2DHG contact always remains grounded during operation. The 2DEG is held at Vil, which accounts 
for the difference in the electron and hole Fermi levels and determines E/^. The Vil ends up being less (~ 35 meV) 
than the GaAs bandgap energy (~ 1.51 eV), due to the other field sources, Vtg, Vbg and the carriers in each well. 
Ideally, n and p are controlled only by their nearest gate, the top and bottom gates, respectively, due to screening. 
However, the system is over-determined (two densities and three voltages) so the same densities can be achieved at 
different gate voltage settings. With the 2DEG held at Vil with respect to ground, all the circuitry connected to it 
must also be held at Vil, necessitating the use of an isolation transformer to break the ground of the signal source. 

The n and p in the uEHBL were set by adjusting Vtg, Vbg and Vil and measured using low- field Hall measurements. 
To characterize transport, the resistivity p was measured in each layer as a function of n and p, E/£ and temperature 
T. The mobility in each layer was calculated from the resistivity and density according to p p = l/pep p and p n = 
l/nep n . Resistivity and Hall measurements were made by standard van der Pauw methods using low- frequency, 
lock-in technique with separate 20 nA excitation currents in both layers. Coulomb drag measurements were used to 
estimate the change in rf as V;l changed. For these measurements a 10 nA current was driven in the 2DEG, while 
the induced voltage in the 2DHG was measured with a high-impedance detection circuit. The constant temperature 
measurements were all taken at T = 0.3 K in a He 3 refridgerator. 

IV. RESULTS & DISCUSSION 

To establish some basis for comparison, the mobility of each layer is plotted in Fig. [5] as a function of its density, 
n or p, at different densities in the adjacent layer with Vil = —1-44 V. The resulting interlayer electric field E/£ 
(ks 96 kV/cm) was below the expected breakdown-field in Al.9Ga.1As, equal to 500 kV/cm for low temperature.— In 
the plots of Fig. [5^ an increase in hole mobility p p with increasing n is visible with a weakening dependence as p 
increases. The plots of p p range from .25 to 3.25 x 10 5 cm 2 /V s for p — 4.0 to 13.0 x 10 10 cm -2 . For m* = .45 m e 
this led to a hole scattering time t p ranging between m 6.4 and 83.2 ps. The p p did not have an ideal power law 
p oc p ap dependence on p across this range. However, the conductivity a p used to calculate p p was fit (not shown) to 
a percolation model (a p oc (p — p c ) lp ). The fit parameter 7 P increased from ss .98 to 1.12 and the critical hole density 
p c decreased from 4.1 to 3.0 x 10 10 cm~ 2 as n increased from 2.0 to 19.0 x 10 10 cm~ 2 . This percolation model was 
previously used to describe the 2D metal-insulator transition (MIT) in undoped heterostructures^ This range of p c 
values are below all the p values used in this study so a 2D MIT was not expected in the 2DHG. 

In contrast, in the plots in Fig. f2]o the electron mobility p n shows no dependence on p, the density in its adjacent 
well. The p n ranges from 3.5 to 10.25 x 10 5 cm 2 /V s for n = 3.0 to 13.0 x 10 10 cm" 2 . For m n = .067 m e , the p n led 
to t„ between ps 13.0 and 39.0 ps. In comparison with p p , p n was oc n Q ™ with a n — 1.02; this fit used a zero-density 
offset po = 1-9 x 10 5 cm 2 /V s. Fitting the conductivity a n (not shown) to a similar percolation transport model as 
above, a fit parameter j n sa 1.8 was found for all p. The a n fitting, however, indicated the 2DEG will never enter an 
insulating regime (n -C n c ). 

In Fig. [3J the mobility layer interdependence is more closely investigated. Fig. [3Jl shows a monotonic relationship 
existed between p p and n, the density in the adjacent well. Furthermore, this dependence apparently becomes weaker, 
illustrated by the visible decrease in slope between the datasets, as p increased from 5.0 to 10.0 x 10 10 cm -2 . This 
decrease in slope correctly corresponds to the data in Fig. [2^, where the spread between plots was much larger for 
smaller p. The inset plot of Fig. [3Ji shows the change in backgate voltage AVbg required to maintain a constant 
p while n was increased from 3.0 to 13.0 x 10 10 cm -2 using Vtg- The traces in Fig. (3Jo show that p n was roughly 
independent of p, confirming what was apparent in Fig. ^Bp. Any non-linearity in this data was attributed to noise in 
the measurements. 

A forthcoming, full calculation for mobility-density data on this structure was qualitatively consistent with transport 
in each layer being dominated by a uniform background (channel) impurity density^, which the results above are 
also consistent with. The increasing p p in Fig. [3Jl would result from increased squeezing of the hole wavefunction. 
Squeezing is known to reduce channel impurity scattering in undoped heterostructures in this density range (< 10 11 
cm" 2 )3 In the uEHBL, the hole wavefunction squeezing increases with increasing n because Vbg, which tilts the 
2DHG confinement potential, was simultaneously decreased (made more negative) to maintain a constant p, as shown 
in the inset of Fig. [3Ji. Furthermore, the decrease in steepness of the slopes of each dataset in Fig. [3Jl as p was 
increased is also consistent since the change in p p is proportional to —AVbg- 

If channel impurity scattering also limited the 2DEG mobility then the results in Fig. [3Jo, that p n is roughly 
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FIG. 2: Mobilities (a) \i v and (b) [in as a function of p and n, respectively, and the adjacent well density at T — 0.3 K and 
V IL = -1.44 V. 



independent of p, imply squeezing of the electron wavefunction was weaker than squeezing of the hole wavefunction 
as density in the respective, adjacent well increased. A direct comparison of the changes in the top and bottom gate 
voltages is to a large degree voided by the device's asymmetry with regard to the cladding layer widths and relatively 
different gate leakages, which was a function of the type of gate and contact combinations for either 2D system (see 
Fig. QJi). However, the mobility data in Fig. [2]still points to channel impurity scattering limiting the 2DEG transport. 
From the power-law dependence of (j. n on n, the a n = 1.02 is clearly below the a = | expected if remote impurity 
scattering is limiting the 2DEG. 1 ? Furthermore, according to the data in Fig. [2] the 2DHG scattering time r p grows 
larger than the 2DEG scattering time r n above n = p>5x 10 10 cm~ 2 . For example, a mobility ratio \i n l li p ~ 3.3 
exists at n = p = 7.0 x 10 10 cm~ 2 . This leads to a scattering time ratio t p /t„ k 2.1. If channel impurity scattering 
limits fj, n than this ratio suggests the 2DEG wavefunction is much wider than the hole wavefunction (see eqns. 4-7 
ir>i£) since that would increase the differential cross section for scattering. A wider 2DEG wavefuncion is expected 
since m* < m*. 
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FIG. 3: Mobilities (a) /x p and (b) fi„ as a function of adjacent well carrier density, n and p, respectively, at T = .3 K and 
Vil — —1.44 V. The inset plot shows AVbg as a function of p for the data in (a). 



In Fig. Q] the same /x p data at p = 5 x 10 10 cm~ 2 and Vil = —1-44 V from Fig. is plotted alongside similar 
measurements at Vil = —1-45 V and —1.43 V. Measurements of \x n under similar conditions (not shown) were also 
taken, but showed no discernable dependence on Vil- The inset plot shows the AVbg required to maintain constant 
p while increasing n from 3.0 to 12.0 x 10 10 cm" 2 at each Vil- With the 2DHG held at ground, making Vil, the 
voltage dropped across the barrier, less negative pulled the 2DEG energy level down, thereby increasing En and 
also ostensibly decreasing interlayer separation d (see Fig. [S]). Based on the prior discussion above, the slopes of the 
datasets in Fig. [4] suggest the largest AVbg would have occurred at Vil = —1.43 V since it has the steepest slope 
and that AVbg increases with decreasing Vil- Both likelihoods are confirmed by the data in the inset plot. 

Following the current argument, the increase in /i p as Vil decreases for fixed n in Fig. 3] also implies increased hole 
wavefunction deformation. However, the reverse, an increase in Vbg &s Vil decreased, was observed. For example, 
at n = 3 x 10 10 cm" 2 a Vbg = -1.647 V resulted for V IL = -1-43 V, which was less than V BG = -1.5735 V at 
Vil — —1.45 V. This result suggests the effect of Vil on squeezing the hole wavefunction overcomes the reverse effect 
from an increasing Vbg to enhance /i p and that narrowing the analysis to only Vbg's influences is insufficient to fully 
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FIG. 4: Mobility fi p as a function of n for Vil = —1.43 V to —1.45 V and (inset) AVbg as a function of Vil at T = .3 K. 
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FIG. 5: The hole drag resistivity pdrag as a function of T at matched density n — p ~ 5.0 x 10 cm 2 for various Vil and 
(inset) schematic of the drag measurement. 



account for the device's behavior. 

The 2DHG Coulomb drag p ( irag measurements, shown in Fig. [5l were taken as function of T at various Vil f° r 
matched density n = p ~ 5.0 x 10 10 cm -2 . At T — 0.3 K a Pdrag ~ -1 ^/sq was equivalently measured for each 
Vil- Using pdrag = rn*/e 2 pTh^e the time it takes for a hole to transfer its momentum to an electron is Th^ e ~ 313 
nsi22 This was much longer than the hole scattering time r p , which varies from ~ 12.7 to 33.2 ps as Vil decreases at 
n = p = 5 x 10 10 cm -2 in Fig. [4] 

The pdrag results in Fig. O also provide some indication that decreasing Vj l squeezes the hole wavefunction and, 
thus, how despite the increase in Vbg the fx p still increased with decreasing Vil at fixed n in Fig. 2J Above T = 1.5 K, 
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FIG. 6: The (a) p v and (b) p e as a function of T at various p and n, respectively, for Vil = —1.44 V. 



Pdrag decreased as Vil was decreased, which is expected, based on theory, if d increased^ To analyze this result the 
ratio (pdrag(A)/ pdrag(B)) 1 /4 was determined, where A and B were the various Vil and A > B. For constant density 
this ratio is proportional to the ratio of interlayer separation d(B) /d(A) at each Vi l ■ From the ratio calculation, the 
decrease AVu = — 10 mV led to an increase Ad « 5%. For a nominal separation d — 38 nm this equates to 1.9 nm. 
This increase in d is suggestive of the hole wavefunction moving away from the 2DEG and closer to the edge of the 
confinement potential causing squeezing to occur. Thus, a decrease in Vil also suggests an increase in squeezing and 
a commensurate increase in fi p . 

Temperature sweeps of p p (T) and p n {T) as a function of n and p were also taken. The p P (T) results, in Fig. [6^,, 
demonstrate qualitative changes towards an insulator state as p decreased, but not a full transition. Recall, p c ranging 
from « 3.0 to 4.0 x 10 10 cm~ 2 was expected from percolation model fittings of <r p versus p. The p n (T) results in 
Fig. show qualitatively more metallic behavior and the disappearance of the peak at T = 4.2 K as n increased. 
Corresponding temperature sweeps of p n {T) and p P (T) as a function of the adjacent well density, p and n, respectively, 
were also done (not shown). For T < 4 K, the former reflected the \i n data at T = 0.3 K in Fig. [5Jd. For T > 4 K, 
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however, p n (T,p) became smaller as p increased. The p p (T,n) displayed no qualitative changes as n increased from 
4.0 to 12.0 x 10 10 cm -2 aside from a small increase in amplitude, which reflects the increase of p p as n increased in 
Fig. Hi. 



V. CONCLUSION 



The layer interdependence of low temperature transport in a 30 nm barrier uEHBL device was investigated. An 
increase in fi p with increasing n was observed at various p, while no change in p n with increasing p was apparent at 
any n. The former appears to be only an indirect effect, however, since Vbg was simultaneously decreased to maintain 
constant p while n was increased. This led to further tilting of the 2DHG confinement potential, thereby increasing the 
hole wavefunction squeezing and, commensurately, reducing the dominant background impurity scattering. Decreasing 
Vj£ was observed to increase both p p and d. A Ad as 5% increase with AVil = — 10 mV was determined by pdragiT) 
measurements. Finally, measurements of p P {T) and p n {T) showed qualitative hints of a 2DHG transition to an 
insulator state that disappeared as p increased and of the 2DEG becoming more metallic as n increased, respectively. 
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